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Abstract 

We compare and analyse Bose-Einstein correlation data for Tr+p and K'^p collisi- 

m 

ons at ^/s = 22 GeV obtained by the NA22-Collaboration and data for pp collisions at 
\l ^/s = 630 GeV obtained by the UAl-Minimum-Bias-Collaboration. Using a parame- 

trization for a longitudinally expanding source, we observe that correlation lengths and 
radii extracted from fits to the data change significantly if one goes from the NA22- 
to the UAl-data: the transverse radius of the chaotic source increases by about 50% 
while the correlation lengths in longitudinal and transverse directions decrease by 
about 50% and 30%, respectively. The chaoticity parameter remains approximately 
constant or increases. 
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1 Introduction 



Particle interferometry constitutes an important method to obtain information about 
space-time aspects of multiparticle production, as well as information about the quantum 
statistical nature of the source (coherence versus chaoticity), in high energy collisions. 
Since the majority of secondaries produced in such reactions are mesons, one usually con- 
siders Bose- Einstein correlation (BEC) data. In particular, one may study the dependence 
of radii, lifetimes, correlation lengths and the chaoticity of the particle source on the type 
of reaction and on the collision energy ^/s. This is important, among other things, for 
the determination of energy densities in the search for the quark-gluon-plasma. The re- 
sults of such a comprehensive analysis of BEC data may provide useful constraints for 
phenomenological models of multiparticle production. 

Note that in practice the extraction of space-time parameters and chaoticity from 
experimental data is not at all straightforward. To see why that is so, consider, for 
instance, the two-particle BEC function 



where pi(k) and /02(^i7^2) are the one- and two-particle inclusive spectra, respectively. In 



particles, i.e., on 6 independent variables (5 variables for a central collision). Because of 
limited statistics, the experimental data are often analysed and presented as a function of 
the single variable = —q^q^, where q^ = — ki^ is the difference of the four-momenta 
of the particles. Compared to other single variables (such as the rapidity difference or 
the difference of transverse momenta) has two important advantages: (a) it is Lorentz 
invariant and (b) for = 0, the two particles have equal momenta and therefore the 
BEC is maximal. The representation of the BEC data as a function of one variable only 
implies an integration over the remaining momentum components. In the variable Q^, the 





the most general case, C2{ki,k2) depends on the three-momenta ki (i = 1,2) of the two 
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two-particle BEC function reads 

CMQ') = -j^ f (2) 

dur du2 Pi{ki)pi{k2) S [Q^ + (fcf - k'^ f] 

where dijj = d^k / {2'k)'^2E is the invariant phase space volume element. 

Thus, in order to obtain information about the space-time characteristics and chaoticity 
of the source from an analysis of the BEC data as a function of Q^, it is necessary to 
explicitly perform the phase space integrations in (2) and to carefully take into account 
the kinematical cuts of the experiment under consideration^. In Refs. [1, 2], this procedure 
was applied to the new BEC-data in pp collisions at ^/s = 630 GeV obtained by the UAl- 
CoUaboration [3], which cover a range in down to values as low as 10""^ GeV^ . 

Recently, the NA22- Collaboration has measured the two-particle correlation in the 

variable Q'^ for Q'^ > 10"^ GeV'^ for vr+p and K+p collisions at = 22 GeV [4]. Since 

both experiments cover the same wide range in Q, it becomes possible for the first time 

to compare BEC data for hadronic collisions at different CM energies and to extract 

information about the change in the space-time parameters and chaoticity of the particle 

source. It is the purpose of the present paper to do this. In section 2, we briefiy introduce 

the formalism and our model for a longitudinally expanding source. The fits to the data 

and the corresponding differences in the source parameters obtained for the NA22- and 

the UAl-data are presented in section 3. Finally, we discuss our results and compare them 

to the results of other theoretical and experimental investigations in section 4. 

'^Note that these integrations may considerably change the shape of the correlation function. In Refs. 
[1, 2] it was demonstrated that a Gaussian form of C2{ki, ^2) can lead to a correlation in which looks 
like a power law. 
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2 The Expanding Source Model 

In this section we briefly review the quantum statistical formalism and the model of a 
boost-invariant longitudinally expanding source^ '"^ which were derived in [5, 6] and applied 
to the analysis of UAl-data in [1, 2]. 

It is assumed that the system can be described by a Gaussian density matrix"*. In this 
case, all multiparticle distributions can be expressed in terms of two functions I(k) and 
D(ki,k2). In particular, the one- and two-particle inclusive distributions take the form 

p,(k) = D(k,k) + \I(k)f (3) 

piiM.h) = pi(M)pi(k2) 

+ 2Re Diki,k2)I*iki)Iik2) + Diki,k2)Dik2,ki) (4) 

To describe the space-time properties of the source we shall work in the current formalism 
[9]. The functions I(k) and D(ki,k2) are the on-shell Fourier transforms of the coherent 
component of the current, I(x), and the chaotic current correlator, D(x,x'), respectively, 
which can be written as [5, 6] 

I{x) = Ux) (5) 
D{x,x') = fch(x)C(x-x') U(x') (6) 

Here, fch{x) and fc{x), are the space-time distributions of the chaotic and of the coherent 

component, respectively, and C{x — x') is the primordial correlator which describes the 

correlation for an inflnitely extended source. 

^Such a model is based on the scenario of an inside-outside cascade and the presence of a plateau in 

the central rapidity region. 

^We have also considered a model of a "static" (non-expanding) source [5, 6] and applied it to the new 

BEC data. It turns out that the static source model fails to describe the UAl-data. Our results obtained 

for the static source will be published in [7]. 

*The Gaussian form of the density matrix follows from the central limit theorem in the case of a large 

number of independent source elements[8]. 
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For the model of a longitudinally expanding source [6], we introduce the variables 

1 , Xo + 

- m 

2 xq — x\\ 



;q - Sn, 7] = ;^ln— — - (7) 



where r is the longitudinal proper time and r] the space-time rapidity. For the sake of 
mathematical simplicity and in order to keep the number of parameters as low as possible, 
we consider for the quantities fch{x), fc{x) and C{x — x') forms which are invariant under 
boosts of the coordinate frame in longitudinal direction: 

fch{x) OC ^(r - To) exp 1^-— (8) 

fc{x) oc 8{t - Tq) (9) 



C(r] — r]',xj_ — x'j_) oc exp 



Ll V 2 ; 2Ll 



:io) 



The parameters for the expanding source are^ the proper formation time tq, the trans- 
verse radii of the chaotic and the coherent source, Rch,± and Rc,±, the correlation lengths 
in longitudinal and transverse direction, and Lj_, and the value po of the chaoticity 
parameter at some fixed momentum (e.g., at kj_ = 0, p(kj_ = 0) = po). As was discussed in 
detail in [2], the number of free parameters can be reduced by relating them to the shape of 
the one-particle inclusive distribution pi(k). For the expanding source, the remaining four 
free parameters are[2] tq, Rch,±, and p (where it has been assumed that the transverse 
momentum spectra of the chaotic and the coherent component have the same shape, in 
which case the chaoticity p becomes independent of fc, cf. [1, 2]). 



^We have made the simplifying assumptions that the chaotically and the coherently produced particles 
are emitted at the same proper time To and that there are no fluctuations in To (for the more general case, 
see Ref. [2]). 



5 



3 Fits to the Data 

Fig. la shows the UAl-data^ for the two-particle BEC as a function of together with 
the fits obtained in Refs. [1, 2] for a longitudinally expanding, purely chaotic source 
(p = 1). The corresponding values of tq, Rch,± and i^, and the values of x^/rf.o./., are 
listed in Table 1. 

The NA22-data are plotted in Fig. lb. Clearly, the BEC function C2iQ^) at = 
22 GeV is broader than the one at 630 GeV. Naively, one might conclude from this 
observation that the radii and/or lifetimes of the source increase with CM energy. However, 
one must bear in mind that the correlation functions were obtained in different phase space 
regions i7 which enter into the integrations in (2): for the NA22-data, i7 is defined through 
the cuts \y\ < 2, kj_ > 0, whereas for the UAl-data one has \y\ < 3, kj_ > 0.15 GeV. 
It turns out that, if one uses for the NA22-data the same parameter values as for the 
fit to the UAl-data (taking into account the different phase space cuts), one obtains 
unacceptable x^/rf.o./. values, i.e., one fails to reproduce the NA22-data. This implies that 
the observed difference in the widths of C2(Q^) in the NA22- and in the UAl-experiment 
cannot be explained through the difference in the rapidity and transverse momentum cuts 
alone. Rather, it refiects a genuine difference in the geometry and quantum statistical 
properties of the particle sources. To describe the NA22-data, one needs to choose a 
set of parameter values different from the ones used to fit the UAl-data. The resultant 
space-time parameters are listed in Table 1, and the corresponding correlation function is 
compared to the NA22-data in Fig. lb. 

Up to this point, we have assumed that the particles are emitted from a purely chaotic 

source. Before we can discuss the quality of the fits and the s-dependence of the quantum 
''We have normalized both the UAl-data [3] and the NA22-data [4] to unity at = 1 GeV^ in 
order to subtract the elFects of long range correlations which were shown to play an important role at 
UAl-energies[10]. 
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statistical parameters of the source, we have to address the question if the data allow any 
conclusions concerning the presence of a coherent component and its energy dependence. 
It turns out that for the NA22-data fits with acceptable x^/rf.o./. can be obtained for 
a wide range of chaoticity parameters, 0.3 < p < 1. For the UAl-data, the range is 
0.8 < p < 1. Fortunately, the values of the other fit parameters (tq, Rch,± and i^) do not 
strongly depend on the value of p. To illustrate this, we present the results of the best fit 
to the NA22-data for p = 0.6 in Table 1. 

Of course it would be interesting to compare the BEC data at ^/s = 22 GeV and at 
.yi = 630 GeV in the same momentum space region. As was mentioned above, the present 
kinematical cuts are \y\ < 2, kj_ > for the NA22-data and \y\ < 3, kj_ > 0.15 GeV for 
the UAl-data. We would therefore urge both experimental groups to re-analyze their data 
and present C2(Q^) in the restricted phase space region \y\ < 2, k± > 0.15 GeV. This 
would allow to explicitly see the s-dependence of the correlation function without having 
to take into account the effects of different phase space cuts. 

Note that here we have not treated separately the contributions from the different 
resonance decays''. Thus the radii and lifetimes extracted from the BEC data and listed 
in Table 1 must be regarded as effective radii and lifetimes in so far as they represent 
an average of the direct pion component and the components related to pions produced 
from resonance decays. We have also not attempted to correct for the effects of particle 
misidentification. 

4 Discussion of the Results 

We have analysed the new BEC data obtained for vr+p and li'+p collisions at ^/s = 22 GeV 

by the NA22- Collaboration and the data for pp collisions at ^/s = 630 GeV obtained by 
formalism which explicitly takes into account the effect of resonance decays on BEC is presented in 
Ref. [ff]. 
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the UAl-Minimum-Bias-CoUaboration. It was shown that the difference in the NA22- 
and the UAI-data cannot be accounted for by the different phase space cuts of the two 
experiments. 

As can be seen from Table 1, the fits to the data based on the model of a boost-invariant, 
longitudinally expanding source show that the fit parameters change significantly from the 
NA22- to the UAI-data. The transverse radius of the chaotic source, Rch,±, increases by 
about 50% (from ~ 0.7 fm to 1 fm) and the correlation length in longitudinal direction 
Lfj is reduced by more than 50% (it drops from 0.9 fm to 0.4 fm). On the other hand. 
To is found to be almost unaffected (it changes from 1.9 fm to 2 fm). The correlation 
length in transverse direction Lj_ (which in the model considered here can be expressed 
in terms of Rch,± and the average transverse momentum of particles, cf. [2]) decreases 
from ~ 0.8 fm to 0.6 fm. If there is a coherent component, the transverse radius of the 
coherent source remains constant at a value of Rc,± ~ 0.5 fm. Concerning the chaoticity 
parameter one can conclude that it either increases or remains approximately® constant. 

The increase in the transverse radius Rch,± from the NA22- to the UAI-data may be 
related to the difference in the cross sections for the reactions vr+p/li'+p at ^/s = 22 GeV 
and pp at ^/s = 630 GeV , which is due to the dependence of the cross sections both on 
the reaction type and the collision energy. For the total cross sections, one has [12] 



Thus, if one assumes that the transverse radius scales like the square root of the total 

cross section, one would expect an increase of Rch,± by about 60% which is compatible 

with our results (cf. Table 1). For a hydrodynamically expanding source the s-dependence 

of the longitudinal and transverse fiow components may also contribute to the increase of 

the transverse radius. As was shown in Ref. [13], for reasonable values of the freeze-out 
The data are also consistent with a slight decrease from p ~ 1 at = 22 GeV to p ~ 0.9 at 
^ = 630 GeV. 



aPP(22 GeV) 
a^^P{22 GeV) 



1.7, 



ctPP(630 GeV) 
aPP{22 GeV) 



1.5 



(11) 



8 



temperature transverse flow can lead to a decrease of the transverse source radius extracted 
from BEC data (cf. Fig. 3 in [13]). At higher energies, the longitudinal expansion 
dominates and transverse flow becomes less important (cf. Ref. [14]). Therefore, the 
transverse source radius can be expected to increase with the CM energy. 

The correlation lengths reflect dynamical properties of the hadronic medium rather 
than the geometry of the source and thus should not depend strongly on the reaction 
type. The data thus indicate a signiflcant decrease of the correlation lengths and Lj_ 
with increasing CM energy. This is the flrst time that information about the s-dependence 
of the correlation length has been extracted from BEC data ^. 

On the theoretical side, there exist surprisingly few predictions concerning the depend- 
ence of BEC on the CM energy ^/s. A decrease of the width of the BEC function in 
longitudinal direction, C2iq\\), with increasing s was predicted by Bowler in the context of 
a string model[16]. A similar behaviour is also expected in a hydrodynamic scenario: for 
a longitudinally expanding source, the width of C2iq\\) is approximately inversely propor- 
tional to the freeze-out proper time Tf [17]. If the formation time tq of the thermalized 
matter is assumed to be s-independent, Tf increases with the initial energy density which 
in turn increases with s. Consequently, the width of C2iq\\) is expected to decrease with 
increasing CM energy. 

Note that our results concerning the s-dependence of the space-time parameters of 

the expanding source are in qualitative agreement with the experimental observation 

[18, 19, 20, 21] that at flxed ^/s the width of the BEC function decreases with incre- 
^In Ref. [15], an analysis of the dependence of the normalized factorial moments of the multiplicity 
distribution on the width of the rapidity interval for NA22- and for UAl-data led to the conclusion that 
the correlation length in rapidity increases with s. However, in [15] no distinction between long range and 
short range correlations was made. As was shown recently in [10], the long range fluctuations increase 
strongly with s, both in intensity and range. Thus there is no contradiction between the results of [15] 
and those of the present paper. 
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asing rapidity density dn/dy. Fluctuations in dn/dy correspond to fluctuations in the 
hadronic energy W, i.e., the part of the CM energy that goes into the production of se- 
condary particles. Since the average dn/dy increases with s, the dn / dy-dependence of 
the BEC data suggests that the width of the correlation function decreases, which in the 
expanding source model is reflected in an increase of Rch,± and a decrease of i^, Lj_ with 
increasing CM energy. On the other hand, the data show that the "intercept" of the 
correlation function, as obtained from extrapolating a simple Gaussian form to Q = 0, 
considerably decreases with increasing dn/dy, which could be taken as an indication that 
the chaoticity decreases with W, and hence with ^/s. Such a behaviour (i.e. a signiflcant 
decrease of p with increasing s) is not consistent with the results of our present analysis. 
The resolution of this apparent discrepancy may be linked to the fact that the extraction 
of a reliable value for the chaoticity parameter requires knowledge of the behaviour of the 
BEC function near the origin, i.e., at small momentum differences. It would therefore be 
useful also from this point of view if the experimentalists could measure C2{Q^) at flxed 
dn/dy down to values of Q ~ 30 MeV . This kind of analysis would help to determine the 
dependence of the quantum statistical parameters on the two variables dn/dy and s. 

This work was supported by the Federal Minister of Research and Technology under 
contract 06MR731 and the Deutsche Forschungsgemeinschaft (DFG). B.R. Schlei acknow- 
ledges a DFG fellowship. We are indebted to B. Buschbeck (UAl Collaboration), W. 
Kittel (NA22 Collaboration) and U. Ornik for many instructive discussions. 
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Figure Captions 



1 Two particle correlation function C^iQ^), for (a) K'^p/ir'^p collisions at ^/s = 
22 GeV and (b) pp collisions at ^/s = 630 GeV . The fits to the UAl-data [3] and 
the NA22-data [4] were obtained for a purely chaotic, boost-invariant longitudinally 
expanding source. 
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Table Captions 



Table 1 Fit parameter values and / d.o.f. for the expanding source model. 
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Expanding source 


[GeV] 


P 


To [fm] 


Rch,± [fm] 


Lri [fm] 


xVd.o.f. 


630 


1.00 


9 19+0.07 

^•-■-^-0.08 


J^-^'J-o.os 


0.42i°:°^ 


49.77/43 


22 


1.00 


-I OO+0.07 


63+°-°'* 


0.88i°:°^ 


36.97/36 


22 


0.60 


^•^^-0.13 


r, 70+0.03 

U. ' «J_o.04 


o.9ot°:°^ 


31.76/36 



Table 1 
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Figure 1a 
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Figure 1b 



